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ABSTRACT 

We present new analysis of 1 1 intermediate-age (1-2 Gyr) star clusters in the Large Magellanic Cloud based 
on Hubble Space Telescope imaging data. Seven of the clusters feature main sequence turnoff (MSTO) regions 
that are wider than can be accounted for by a simple stellar population, whereas their red giant branches indicate 
a single value of [Fe/H]. The star clusters cover a range in present-day mass from about 1 x 10 4 M to 2 x 10 5 
M Q . We compare radial distributions of stars in the upper and lower parts of the MSTO region, and calculate 
cluster masses and escape velocities from the present time back to a cluster age of 10 Myr. Our main result 
is that for all clusters in our sample with estimated escape velocities v esc > 15 kms" 1 at an age of 10 Myr, the 
stars in the brightest half of the MSTO region are significantly more centrally concentrated than the stars in 
the faintest half and more massive red giant branch and asymptotic giant branch stars. This is not the case for 
clusters with v esc < 10 kms" 1 at an age of 10 Myr. We argue that the wide MSTO region of such clusters is 
mainly caused by to a ~ 200-500 Myr range in the ages of cluster stars due to extended star formation within 
the cluster from material shed by first-generation stars featuring slow stellar winds. Dilution of this enriched 
material by accretion of ambient interstellar matter is deemed plausible if the spread of [Fe/H] in this ambient 
gas was very small when the second-generation stars were formed in the cluster. 
Subject headings: globular clusters: general — Magellanic Clouds) 



1. INTRODUCTION 

For several decades, the standard paradigm for globular 
clusters (GCs) was that they consist of stars born at the same 
time out of the same material. This scenario has faced seri- 
ous challenges over the last decade. It is now known that the 
most massive GCs in our Galaxy such as lu Cen and M 54 host 
multiple red giant branches (RGBs) due to populations with 
different [Fe/H] (e.g., Sarajedini & Layden 1995; Lee et al. 
1999; Hilker & Richtler 2000; Villanova et al. 2007; Carretta 
et al. 2010). Somewhat less massive Galactic GCs such as 
NGC 2808, NGC 1851 and 47 Tuc show multiple sub-giant 
branches (SGBs) and/or multiple main sequences, which are 
typically interpreted as populations with different Helium 
abundances (e.g., Piotto et al. 2007; Milone et al. 2008; An- 
derson et al. 2009). While lower-mass Galactic GCs typically 
do not show clear evidence for multiple populations from op- 
tical broad-band photometry, spectroscopic surveys do show 
that light elements such as C, N, O, F, and Na show signif- 
icant star-to-star abundance variations (often dubbed "Na-O 
anticorrelations") within all Galactic GCs studied to date in 
sufficient detail (Carretta et al. 2009, and references therein). 
Since these abundance variations have been found among 
RGB stars as well as main sequence (MS) stars within a given 
GC (Gratton et al. 2004), the suggested cause of the varia- 
tions is that secondary generation(s) of stars formed out of 
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material shed by an older, evolved population within the clus- 
ter. While the chemical processes responsible for causing the 
light element abundance variations have largely been identi- 
fied as proton-capture reactions in Hydrogen burning at high 
temperature (> 40 x 10 6 K, see e.g. Gratton et al. 2004), the 
old age of Galactic GCs has precluded a clear picture of the 
time scales and hence the types of stars involved in the chem- 
ical enrichment of the second-generation stars. Currently, the 
most popular candidates are ( i) intermediate-mass AGB stars 
(4 < M /M < 8, hereafter IM-AGB; e.g., D' Antona & Ven- 
tura 2007 and references therein), ( ii) rapidly rotating massive 
stars (often referred to as "FRMS"; e.g., Decressin et al. 2007) 
and (Hi) massive binary stars (de Mink et al. 2009). 

Recently, deep CMDs from images taken with the Ad- 
vanced Camera for Surveys (ACS) aboard the Hubble Space 
Telescope (HST) provided conclusive evidence that sev- 
eral massive intermediate-age star clusters in the Magellanic 
Clouds host extended and/or multiple main sequence turn- 
off (MSTO) regions (Mackey et al. 2008; Glatt et al. 2008; 
Milone et al. 2009; Goudfrooij et al. 2009, hereafter Paper 
I; Goudfrooij et al. 2011, hereafter Paper II), in some cases 
accompanied by composite red clumps (Girardi et al. 2009; 
Rubele et al. 2011). To date, these observed properties have 
been interpreted in three main ways: ( i) Bimodal age distribu- 
tions (Mackey et al. 2008; Milone et al. 2009), (ii) age spreads 
of 200-500 Myr (Paper II; Girardi et al. 2009; Rubele et al. 
2010, 201 1), and (Hi) spreads in rotation velocity among turn- 
off stars (Bastian & de Mink 2009). 

In this third paper of this series we study the dynamical 
properties of 11 intermediate-age star clusters in the LMC: 
seven clusters that exhibit extended MSTO regions (hereafter 
eMSTOs) and four that do not. We determine and compare 
radial distributions of cluster stars at different evolutionary 
phases and evaluate the evolution of the clusters' masses and 
escape velocities from an age of 10 Myr to their current age. 
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This analysis reveals new findings relevant to the assembly of 
these intermediate-age star clusters and their evolutionary as- 
sociation with multiple stellar populations in ancient Galactic 
globular clusters. 

The remainder of this paper is organized as follows. § 2 
presents the cluster sample. In § 3 we determine the radial 
distributions of stars in various different evolutionary phases 
and study the dependencies of these distributions with cluster 
escape velocities as a function of time. § 4 uses our new re- 
sults to constrain the origin of multiple populations in these 
clusters, and § 5 presents our main conclusions. 

2. TARGET CLUSTERS 

Our main sample of intermediate-age clusters is that pre- 
sented in Paper II. All these clusters were observed within one 
HST/ACS program using the same observational setup, using 
both short and long exposures to yield high-quality photome- 
try throughout the CMDs. The sample consists of star clusters 
in the LMC with integrated UBV colors consistent with SWB 
(Searle et al. 1980) parameter values in the range IV B - VI, 
which translates to estimated ages between roughly 1.0 and 
2.5 Gyr. This turned out to be fully consistent with the ages 
actually found from isochrone fitting. The main properties of 
these star clusters, all of which were found to host eMSTO 
regions, are listed in Table 1 . For comparison with clusters 
in the same age range that do not contain eMSTO regions, 
we also use the four LMC clusters NGC 1644, NGC 1652, 
NGC 1795, and IC 2146 that were studied by Milone et al. 
(2009) using HST/ACS photometry. To our knowledge, these 
are the only four clusters in the LMC in the age range 1 - 
2 Gyr that are known to date not to harbor eMSTO regions 
(Milone et al. 2009) from high-quality photometry at HST 
resolution. 

3. RADIAL DISTRIBUTIONS OF STARS IN DIFFERENT REGIONS 
OF THE CMD 

3.1. Motivation 

Since the eMSTO regions in our target clusters may be due 
to the presence of more than one simple stellar population, 
it is important to find out whether different subregions of the 
eMSTO have intrinsically different radial distributions, i.e., 
differences beyond those that can be expected for simple (co- 
eval) stellar populations due to dynamical evolution. A well- 
known example of the latter is mass segregation due to dy- 
namical friction which slows down stars on a time scale that 
is inversely proportional to the mass of the star (e.g., Saslaw 
1985; Spitzer 1987; Meylan & Heggie 1997) so that massive 
stars have a more centrally concentrated distribution over time 
than less massive ones. Another reason for studying the radial 
distributions of stars in relevant evolutionary phases in the 
CMD is that the masses and ages of these clusters are such 
that the two-body relaxation times of stars within the clusters 
are comparable to the cluster ages: The mean two-body re- 
laxation time fj-eiax of stars within a cluster's half-mass radius 
is 

fi-eiax « — — W ~ — — rb 5 (GM c] y - 5 (1) 
81nJV 81nJV h 

(Binney & Tremaine 1987) where f cross is the crossing time 
at the half-mass radius, N is the number of stars in the clus- 
ter, r/, is the half-mass radius, G is the gravitational constant, 
and Ai c i is the cluster mass. f re i ax values for the clusters in 
our sample are listed in Table 2. For a typical cluster among 
the 50% most massive clusters in our sample, N ~ 1.5 x 10 5 , 
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Figure 1. Illustration of the three regions used to derive the radial surface 
number density distributions shown in Fig. 2, superposed onto the CMD of 
NGC 1806 as shown in the middle panel of Fig. 2 of Paper II. The "upper 
MSTO region" is outlined by blue solid lines, the "lower MSTO region" by 
red lines, and the "RC, RGB, and AGB" by black lines. 

Ai c \ ~ 10 5 M Q , and rj, ~ 8 pc, resulting in f le i ax ~ 2 Gyr. 
Hence, dynamical imprints from the formation epoch of the 
stars within these clusters may well be still observable, in con- 
trast to the situation in the vast majority of Galactic globular 
clusters (wCen is the exception, see e.g. Bellini et al. 2009). 
This renders these intermediate-age clusters excellent probes 
of the nature of multiple stellar populations in star clusters. 

3.2. Measurements and Implications 

We derived radial distributions of the following regions on 
the CMDs of the seven star clusters in our sample that contain 
eMSTOs: (i) 'all' stars in the CMD, (ii) the 'upper half of the 
eMSTO region, (Hi) the 'lower half of the eMSTO region, 
and (iv) the Red Clump (RC) stars and the RGB and AGB 
stars above the RC. Note that regions (ii) - (iv) were found 
in Paper II to have less than 5% contamination by field stars 
for every cluster. Fig. 1 depicts these regions on the CMD of 
NGC 1806 with solid lines, while the completeness-corrected 
radial distributions of stars in these regions are shown in Fig. 
2. The size of radial bins in all panels of Fig. 2 was deter- 
mined by the requirement that there be at least 30 stars in the 
outermost radial bin for any of the CMD regions considered. 

Fig. 2 shows several items of interest, discussed in turn be- 
low. 

(i) For many of the brighter (more massive) clusters in our 
sample, the stars in the upper (brighter) half of the 
eMSTO region are significantly more centrally concen- 
trated than the stars in the lower (fainter) half of the eM- 
STO. In fact, the stars in the upper half of the eMSTO 
show a central concentration that is even stronger than 
that of (more massive) stars in the upper RGB/AGB 
in the clusters NGC 1751 and NGC 1846 (and, to a 
somewhat lesser extent, NGC 1806). This means that 
mass segregation cannot be the cause of these differ- 
ences. The radial distribution of the RGB/AGB stars 
in these three massive clusters typically follows that of 
the lower half of the eMSTO except for the innermost 
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Table 1 

Main properties of the eMSTO star clusters studied in this paper. 



Cluster 


V 


Ref. 


SWB 


Age 


[Z/H] 


Av 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


NGC 1751 


11.67 ±0.13 


1 


VI 


1.40±0.10 


-0.3 ±0.1 


0.40 ±0.01 


NGC 1783 


10.39 ± 0.03 


1 


V 


1.70±0.10 


-0.3 ±0.1 


0.02 ± 0.02 


NGC 1806 


11.00 ±0.05 


1 


V 


1.67±0.10 


-0.3 ±0.1 


0.05 ±0.01 


NGC 1846 


10.68 ± 0.20 


1 


VI 


1.75±0.10 


-0.3 ±0.1 


0.08 ±0.01 


NGC 1987 


11.74 ±0.09 


1 


IVB 


1.05 ±0.05 


-0.3 ±0.1 


0.16 ±0.04 


NGC 2108 


12.32 ± 0.04 


2 a 


IVb 


1.00 ±0.05 


-0.3 ±0.1 


0.50 ±0.03 


LW431 


13.67 ± 0.04 


T 


VI 


1.75 ±0.10 


-0.3 ±0.1 


0.15 ±0.01 



Note. — Column (1): Name of star cluster. (2): Integrated V magnitude. (3): Reference of V magnitude. Ref. 1: Goudfrooij et al. (2006); Ref. 2: Bica et al. 
(1996). (4): SWB type from Bica et al. (1996). (5) Age in Gyr from Paper II. (6) Metallicity from paper II. (7) A v from Paper II. 



"uncertainty only includes internal errors associated with measurements of 
cluster and one background aperture. 

few radial bins where it is intermediate between the ra- 
dial distributions of the upper and lower halves of the 
eMSTO region. 

(ii) Interestingly, this difference in central concentration be- 

tween the upper and lower halves of the eMSTO is not 
quite as strong for NGC 1783, the most luminous clus- 
ter in our sample. The case of NGC 1783 is discussed 
in some detail in § 3.3.4 below. 

(iii) The difference in central concentration between the up- 

per and lower halves of the eMSTO is insignificant for 
the lower-luminosity clusters LW431, NGC 1987, and 
NGC 2108. 

(iv) Finally, the radial gradient of 'all' stars in the clusters 

(open circles in Fig. 2) is significantly shallower than 
that of the MSTO and RGB/AGB populations. This 
mainly reflects the relatively strong contamination by 
LMC field stars in the lower MS region of the CMD of 
these clusters (see Paper II). 

We conclude that the upper and lower parts of the eMSTOs 
in (at least some of) these star clusters correspond to intrinsi- 
cally physically different populations which seem to have ex- 
perienced different amounts of violent relaxation during their 
collapse and/or different dynamical evolution effects. 

We believe that this finding is relevant in the context of 
the two main interpretations of the eMSTO regions in the 
recent literature: (i) bimodal age distributions or spreads in 
age (Mackey et al. 2008; Milone et al. 2009; Girardi et al. 
2009; Paper I; Paper II) and (ii) spreads in rotation velocity 
among turnoff stars (Bastian & de Mink 2009, but see Gi- 
rardi et al. 2011). Specifically, if the width of the eMSTO 
in intermediate-age star clusters is mainly caused by one or 
more secondary generation(s) of stars having formed within 
the cluster from material shed in slow winds of stars of the 
first generation and retained by the cluster, the simulations of 
D'Ercole et al. (2008) show that the younger generation (i.e., 
the upper half of the eMSTO region) would be more centrally 
concentrated than the older one if the age of the cluster is 
less than or similar to the half-mass relaxation time (which is 
the case for the four most massive clusters in our sample, cf. 
Table 2). Conversely, if the width of the eMSTO is mainly 
due to a range in rotation speeds among stars with masses 
1.2 < M/Mq < 1.7 (Bastian & de Mink 2009), one might 
expect more rapid rotators (i.e., the lower half of the eMSTO 
region according to Bastian & de Mink 2009) to be initially 



more centrally concentrated than less rapid rotators since ob- 
servational studies have found young stars in dense star clus- 
ters to show higher rotation rates than similar stars in the field 
and hence born in presumably less dense stellar aggregates 
(e.g., Keller 2004; Strom et al. 2005). While detailed mod- 
elling of the dependencies of the radial distribution of stars in 
star clusters on stellar rotation velocity and cluster age is still 
lacking, the observation that the upper (i.e., brighter) half of 
the eMSTO population is significantly more centrally concen- 
trated than the lower (i.e., fainter) half in several star clusters 
in our sample seems to indicate that age effects are responsi- 
ble for the broadening of the eMSTO (at least in those clus- 
ters). 

3.3. Trends with Cluster Escape Velocity as a Function of 

Time 

In a scenario where secondary generations form from mate- 
rial shed in slow stellar winds of the first generation, expected 
wind speeds for intermediate-mass AGB (IM-AGB) stars are 
about 10-20 kms" 1 (e.g., Vassiliadis & Wood 1993; Mar- 
shall et al. 2004). For fast-rotating massive stars (FRMS), 
wind speeds range between about ten to a few hundreds of 
kms" 1 (Porter 1996; Porter & Rivinius 2003; Wiinsch et al. 

2008) . For massive interacting binary stars, one can expect 
wind velocities of a few tens of kms" 1 as well: Observations 
of the well-studied system RY Scuti have shown ejection ve- 
locities of 30-70 kms" 1 (Smith et al. 2002; de Mink et al. 

2009) . If any (or all) of these types of stars provide the mate- 
rial for secondary generations in star clusters, one would ex- 
pect the ability of star clusters to retain this material to scale 
with their escape velocities at the time such stars are present 
in the cluster (i.e., at ages of ~5-30 Myr for the massive 
stars and ~ 50 -200 Myr for IM-AGB stars). Furthermore, 
a well-defined relation has been suggested to exist between 
the most massive star in a cluster and the initial cluster cloud 
mass in that initially more massive clusters host more mas- 
sive stars. The physical effect that causes this relation may 
be that retention of feedback energy from more massive stars 
requires more massive gas clouds (Weidner & Kroupa 2006). 
Hence, if massive stars are a significant source of the material 
used to form secondary generations in clusters, one would ex- 
pect the existence of a relation between the fraction of second- 
generation stars and the mass of the cluster. 

With this in mind, we estimate masses and escape velocities 
of the clusters in our sample for both the present time and at 
an age of 10 7 yr (i.e., after the cluster has survived the era 
of gas expulsion and violent relaxation, see e.g. Baumgardt 
et al. 2008) as follows. For comparison with clusters in the 
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Figure 2. Radial surface number density distributions of stars in different evolutionary phases, as illustrated in the CMD in Figure 1 with solid lines. Cluster 
names are given in the top right in each panel. Open circles: all stars in the CMD. Black circles: Stars in the RC, RGB, and AGB. Red circles: Lower MSTO 
region. Blue circles: Upper MSTO region. The absolute surface number density values on the Y axis refer to the open circles ("all stars"). The surface number 
densities of the other star types are normalized to that of "all stars" at the outermost radial bin. See discussion in § 3. 



same age range that do not contain eMSTO regions, we also 
do these calculations for the four LMC clusters NGC 1644, 
NGC 1652, NGC 1795, and IC 2146 from the study of Milone 
et al. (2009). 

3.3.1. Present-day Masses and Half-Mass Radii 

Present-day masses and half-mass radii for the clusters in 
our sample are adopted from Paper II (see its Tables 1 and 3), 
using Ai /Ly ratios of the SSP models of Bruzual & Chariot 
(2003). For the four clusters found by Milone et al. (2009) 
not to contain eMSTOs, we use the total V-band magnitudes 
in Bica et al. (1996) along with age, distance, and foreground 
reddening values from Milone et al. (2009). Since we cannot 
find half-light radii published for those four clusters in the lit- 



erature, we assume an average value of clusters in the LMC 
within the age range 1-2 Gyr from the literature (Mackey 
& Gilmore 2003; Paper II) along with a suitably large uncer- 
tainty, namely rj, = 8 ± 4 pc. 

3.3.2. Early Cluster Mass Loss 

We make a distinction between model clusters with and 
without initial mass segregation, the inclusion of which can 
have a strong impact on cluster mass loss and dissolution due 
to the strong expansion of the cluster in response to rapid 
mass loss associated with supernova (SN) type II explosions 
of first-generation stars (e.g., Chernoff & Weinberg 1990; 
Fukushige & Heggie 1995; Vesperini et al. 2009). To prop- 
erly represent the case of clusters with initial mass segregation 
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for this paper, a selection needs to be made as to the initial 
properties of the cluster since the recent literature shows that 
different initial properties can lead to a wide range of of clus- 
ter dissolution times (Baumgardt & Kroupa 2007; Baumgardt 
et al. 2008; D'Ercole et al. 2008; Vesperini et al. 2009). We 
consider the evolution of cluster mass and half-mass radius 
in the simulation called SG-C30 in D'Ercole et al. (2008; 
Vesperini 201 1), which involves a tidally limited model clus- 
ter that features a moderate degree of initial mass segrega- 
tion. This simulation is selected for the following reasons: 
(i) mass-segregated clusters of this type can survive the rapid 
SN-driven mass loss era (Vesperini et al. 2009); (ii) it yields a 
ratio of first-to-second-generation stars of ^ 1 : 1 at an age of 
2 Gyr (see Fig. 16 of D'Ercole et al. 2008), similar to (though 
somewhat lower than) what has been found for the more mas- 
sive clusters in our sample (see Paper II; Milone et al. 2009). 
Note that simulations of clusters with initially more strongly 
mass-segregated stellar distributions can yield results similar 
to that of the SG-C30 simulation if the cluster does not fully 
fill its Roche lobe (Vesperini et al. 2009). We also account for 
the effect of mass loss due to stellar evolution, most of which 
takes place in the first 10 8 yr. The evolution of this mass loss 
is derived using the Bruzual & Chariot (2003) SSP models 
using the Chabrier (2003) IMF 5 . The fast stellar evolution- 
driven mass loss in the first 10 8 yr, in particular that due to 
feedback and supernovae from massive O and B-type stars, 
also causes an expansion of the cluster. This is accounted for 
by following the prescriptions of Hills (1980) assuming adi- 
abatic expansion, since the crossing time of star clusters is 
much shorter than the time scale for stellar evolution-driven 
mass loss (see also Baumgardt et al. 2008). 

3.3.3. Long-term Cluster Mass Loss 

At ages > 100 Myr we account for long-term dynami- 
cal cluster destruction mechanisms and follow Fall & Zhang 
(2001) and McLaughlin & Fall (2008) who find that the 
change in shape of the star cluster mass function from young 
(< 5 x 10 8 yr) to old (> 10 10 yr) systems is mainly due to 
mass-density-dependent evaporation (also often called "two- 
body relaxation"). We use the rate of mass loss due to evapo- 
ration derived by McLaughlin & Fall (2008, their Eq. (5)): 



7J^~ 1100 



Ph 



M pc- 



1/2 



M Gyr" 1 



(2) 



where TJeV is the mass loss rate due to evaporation averaged 
over a cluster lifetime and = 3 M c \/%nr\ is the half-mass 
density of the star cluster. The evolution of cluster masses 
is then calculated for an appropriately large grid of masses 
and half-mass radii at an age of 13 Gyr, both with and with- 
out initial mass segregation (implemented as described in the 
previous paragraph). Cluster mass as a function of time for 
individual star clusters in our sample are evaluated by means 
of linear interpolation between grid lines using masses at the 
mean age of the cluster in question. The process of calcu- 
lating cluster masses as a function of time is illustrated in 
Fig. 3 for the case of NGC 1846. For a typical cluster in our 
sample, the difference in initial cluster mass between mass- 
segregated and non-mass-segregated star clusters is of order 
A log A4 C \ w 0.1-0.15 dex, with the initially mass-segregated 
clusters losing a larger fraction of their initial mass. 

5 Use of the Salpeter (1955) IMF results in initial cluster masses that are 
~ 9% lower than those calculated here. 
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Figure 3. Illustration of the method used to calculate cluster masses as func- 
tion of age. (Black) solid lines and the short-dashed line represent calcula- 
tions for initially non-mass-segregated clusters, whereas (red) long-dashed 
lines and the (red) dashed-dotted line represent calculations for an initially 
mass-segregated cluster model (see discussion in § 3.3). Four model calcu- 
lations are shown for the current half-mass radius of NGC 1846 and spaced 
by log(_A/f c i/MQ) = 0.1 at an age of 13 Gyr. The short-dashed and dashed- 
dotted lines depict the adopted mass evolution for NGC 1846. The (blue) 
dotted lines indicate the current age and mass of NGC 1846. 

3.3.4. Escape Velocities 

Escape velocities required to reach the tidal radius of the 
cluster are calculated as follows: 



c(0=/c V^T kms 



(3) 



where A4 c \(t) is the cluster mass in M at time t, r^t) is the 
cluster's half-light radius in pc at time t, and f c is a coeffi- 
cient that takes the dependence of v esc on the concentration 
index c of King (1962) models into account. Values for/ c are 
taken from Georgiev et al. (2009). For convenience, we define 
Vesc.7 = i'esc(10 7 yr) hereinafter. 

To parametrize the relative central concentrations of the up- 
per and lower half of the eMSTO region in the star clusters, 
we define a "relative concentration parameter" as follows: 



'core, low 



^ (^core.upp ^bck.upp) / 

Crei = ■ / 

^bck.upp / ^bck,low 



^bck.low) 



(4) 



where n r 



and «bckupp are surface number densities of 



stars in the upper half of the eMSTO region within the King 
core radius and the background region, respectively, and 
«core,iow and «bck,iow are the same parameters for stars in the 
lower half of the eMSTO region. The background region was 
defined as the outermost radial bin used for the plots in Fig- 
ure 2. Values for A4 C \, r^, t re \. dx , v esc , and C re i are given in 
Table 2, both for the present time and at an age of 10 7 yr. We 
adopt an uncertainty of 20% for present-day cluster mass val- 
ues and use that in calculating uncertainties of f re i a x and v eS c. 

C re i is plotted versus cluster v esc at the current age and at an 
age of 10 7 yr in the two panels of Figure 4. By definition, we 
set C re i = 1 for the clusters found by Milone et al. (2009) not 
to host eMSTOs. The left panel shows that the clusters with 
the highest current escape velocities tend to have the high- 
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Figure 4. Left panel: Parameter C re i, a measure of the relative central concentration of the upper versus the lower half of the MSTO region (cf. Eq. 4), plotted 
against current escape velocity. Right panel: C re i against escape velocity at an age of 10 7 yr. Filled circles represent data for clusters in our sample using a model 
without initial mass segregation while open (red) circles do so using a model involving initial mass segregation. Clusters shown by Milone et al. (2009) not to 
exhibit eMSTO regions are represented by filled (grey) squares and open (green) squares for models without and with initial mass segregation, respectively. The 
names of the star clusters in our sample are indicated to the upper or lower right of their data points. See discussion in § 3.3. 



est value of C re i, although the current escape velocities are 
barely high enough to retain material from even the slowest 
reported winds of IM-AGB stars (i.e., ~ 10 kms" 1 ). How- 
ever, the escape velocities were significantly higher at an age 
of 10 7 yr as illustrated by the right panel of Figure 4. More- 
over, the clusters with v esC 7 > 15 kms" 1 have significantly 
higher values of C le i than those with v esc .7 < 10 kms" 1 , and 
the correlation between C le i and v esc becomes stronger than in 
the left panel. These trends are present whether or not one 
takes initial mass segregation of stars into account in the cal- 
culations of v'esc.7- Summarizing these results, it appears that 
intermediate-age star clusters in the LMC with higher (initial) 
escape velocities have upper ("younger") MSTO stars that 
are more centrally concentrated than their counterparts in the 
lower ("older") MSTO. This finding is discussed further in § 4 
below. The one cluster that seems to deviate slightly from the 
trends mentioned above and shown in Figure 4 is NGC 1783. 
This is the cluster with the the largest radius among our sam- 
ple. In fact, its radial extent is so large that the field of view of 
the HST/ACS image is too small to allow a robust determina- 
tion of its King concentration parameter (see Paper II), which 
translates into a relatively large uncertainty of its half-light ra- 
dius and hence of its escape velocity. Future determinations of 
its rh from wider-field imaging should help improve the accu- 
racy of its escape velocity and hence its position in Figure 4. 

4. IMPLICATIONS REGARDING THE NATURE OF MULTIPLE 
POPULATIONS IN STAR CLUSTERS 

The results described above have implications regarding the 
nature of multiple stellar generations in star clusters in gen- 
eral, including the situation seen in many (ancient) globular 
clusters in our Galaxy. To constrain the scope of this broad 
topic somewhat, we restrict the following discussion to the 
case of star clusters that were never massive enough to retain 
gas expelled by energetic supernova (SN) explosions and/or 
to capture significant numbers of field stars (e.g., from their 
host dwarf galaxies). In practice this restriction corresponds 



roughly to the exclusion of ancient star clusters with current 
masses > 2 x 10 6 M© (cf. Bastian & Goodwin 2006; Fell- 
hauer et al. 2006) at an age of 13 Gyr. The LMC star clusters 
studied in this paper are and were indeed always less massive 
than that limit. 

Star clusters in our sample with initial escape velocities 
Vesc,7 > 15 kms" 1 have C re i values that are significantly higher 
than unity and correlate with v esCj 7. Hence we postulate that 
the extended morphology of their MSTO is mainly caused by 
a range in age, although we cannot formally rule out that the 
widening of the MSTO is partly due to an additional physical 
process for those clusters as well. 

As to the star clusters in our sample with initial escape 
velocities v esC 7 < 10 kms" 1 , they have C re i values consis- 
tent with unity. Under the assumption that the material 
used to form secondary generations of stars within clusters is 
mainly produced internally from stellar winds with v > 10- 15 
kms" 1 , it seems plausible that eMSTO regions in such clus- 
ters may be due in part to physical processes unrelated to a 
range in age, e.g., one or more processes that are independent 
of (or only slightly dependent on) the location of stars within 
the cluster. This could be a range of stellar rotation veloc- 
ities as proposed by Bastian & de Mink (2009) and/or other 
(yet unidentified) physical processes, possibly related to inter- 
actions between close binary stars since the binary fractions 
are significant (15-35%, cf. Paper II) in the clusters studied 
here. Note however that the fact that the clusters found by 
Milone et al. (2009) not to harbor eMSTOs are all at the low 
end of the mass (and v es c,7) range of the clusters considered 
here seems inconsistent with the stellar rotation scenario in 
that one would not a priori expect the distribution of stellar 
rotation rates to depend on cluster mass or escape velocity. In 
this context we note that an age range would still be able to ac- 
count for a widening of the MSTO in low-mass clusters in our 
sample in case their escape velocities at an age of 10 7 yr turn 
out to be higher than our estimates. This would, for example, 
be the case if the initial mass segregation of such clusters was 
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stronger than in the model we used in § 3.3.2. In that case, 
the low values of C re i in these low-mass clusters can still be 
explained by their short half-mass relaxation times which are 
significantly shorter than their age (cf. Table 2) so that initial 
differences in radial distribution between generations of stars 
would now have been eliminated. 

These results impact our understanding of multiple popu- 
lations in ancient globular clusters. For example, the find- 
ing of multiple populations and a Na-O anticorrelation in M 4 
(Marino et al. 2008, and references therein) was described as 
surprising given its low mass of ~ 6 x 10 4 M©. However, 
using its current values for half-light radius = 2.23 pc and 
King concentration parameter c = 39 (Harris 1996), an age of 
13 Gyr, and the methodology outlined in § 3.3, we calculate a 
present-day escape velocity v esc =16 kms -1 and a v esc .7 = 52 
kms -1 for M4 (assuming no initial mass segregation). This 
is more than sufficient to have retained chemically enriched 
material from slow stellar winds of the first stellar genera- 
tion. In fact, for typical half-mass cluster radii < 3 pc, we 
would expect ancient star clusters with current masses as low 
as ~ 5 x 10 3 Mq to have been able to retain such material (see 
also Conroy 201 1 for a scenario based on evolution of cluster 
masses). 

Similarly, recent results on the presence or absence of sig- 
nificant light-element abundance variations within LMC clus- 
ters can be understood by considering the clusters' initial es- 
cape velocities. Mucciarelli et al. (2009) measured Na and O 
abundances for RGB stars in three old metal-poor GCs in the 
LMC which have current masses in the range 2 - 4 x 10 5 M Q 
(see Mackey & Gilmore 2003), and found clear evidence for 
Na-O anticorrelations similar to those observed in Galactic 
GCs. On the other hand, Mucciarelli et al. (2008) measured 
element abundances in four relatively massive intermediate- 
age (1-2 Gyr) clusters in the LMC and found only small vari- 
ations of element abundance ratios within each cluster. How- 
ever, focusing on their [Na/Fe] measurements, we notice that 
two clusters in their sample (NGC 1978 and NGC 2173) ex- 
hibit a range of [Na/Fe] values A[Na/Fe] > 0.4 dex which 
is a > 4cr effect where a is the typical measurement error of 
[Na/Fe] 6 . More recently, Mucciarelli et al. (2011) measured 
element abundances of 14 stars in the young globular cluster 
NGC 1866 in the LMC and found no sign of a Na-O anticorre- 
lation (nor a significant spread in [Na/Fe]). Since they also did 
not detect Na-O anticorrelations in the massive intermediate- 
age clusters in their 2008 paper, they argued for "a different 
formation/evolution scenario for the LMC massive clusters 
younger than ~ 3 Gyr with respect to the old ones". To eval- 
uate whether the presence or absence of significant spreads 
in [Na/Fe] in the LMC clusters studied by Mucciarelli et al. 
(2008; 2009; 201 1) could instead "simply" be due to the clus- 
ters' masses and escape velocities at an age of 10 Myr, we first 
adopt these clusters' current masses, radii, and ages from the 
compilation of Mackey & Gilmore (2003). Since NGC 1978 
was not included in that compilation, we adopt its age, [Fe/H] 
and E(B-V) from Milone et al. (2009), we assume a radius 
of rh = 8 ±4 pc (cf. § 3.3), and we estimate its current mass 
from its integrated V magnitude in Goudfrooij et al. (2006) 
and the M /L v ratios of the Bruzual & Chariot (2003) SSP 

6 The negligible variation of [O/Fe] in these clusters found by Mucciarelli 
et al. (2008) can be explained by their moderately high metallicities for which 
depletion of O in IM-AGB stars is much smaller than for the much more 
metal-poor ancient GCs in the LMC (Ventura & D'Antona 2009; Conroy 
2011). 
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Figure 5. Top panel: A ([Na/Fe]), defined as the measured spread in [Na/Fe] 
among stars in LMC clusters studied by Mucciarelli et al. (2008, 2009, 201 1), 
versus cluster escape velocity at an age of 10 Myr. Error bars for A ([Na/Fe]) 
indicate the typical uncertainty of individual [Na/Fe] measurements. Bottom 
panel: Present-day cluster mass against escape velocity at an age of 10 7 yr 
for the same clusters. Filled (black) circles represent data for clusters with 
A([Na/Fe]) > 0.4 (a 4cr effect) using a model without initial mass segre- 
gation while open (red) circles do so using a model involving initial mass 
segregation. Clusters with A([Na/Fe]) < 0.4 are represented by filled (grey) 
squares and open (green) squares for models without and with initial mass 
segregation, respectively. See discussion in § 4. 

models using the Chabrier IMF. We then estimate these clus- 
ters' masses and escape velocities at an age of 10 Myr us- 
ing the methodology described in § 3.3. We plot the observed 
spreads in [Na/Fe] (defined here as A( [Na/Fe])) from Muc- 
ciarelli et al. (2008; 2009; 201 1) as well as the current cluster 
masses versus v esc .7 in Figure 5. Note that (i) the distinction 
between clusters with and without significant (~4cr) spreads 
in [Na/Fe] is found at v esc .7 w 15 kms" 1 at an age of 10 Myr, 
and (ii) A([Na/Fe]) correlates with v esc 7, indicating that clus- 
ters with higher initial escape velocities are able to retain ma- 
terial with higher variance of light-element abundances. This 
is consistent with the results shown in Figure 4 for the C le i 
parameter of the clusters in our sample, and underlines the 
importance of the clusters' initial escape velocities in their 
ability to retain slow wind material from stars of the first gen- 
eration. 

Our results for the star clusters with v esC) 7 > 15 kms -1 
are generally consistent with the conclusions drawn for 
NGC 1846 in Paper I (see also Renzini 2008; Conroy & 
Spergel 201 1). In particular, we believe our results are most 
consistent with the "in situ" scenario (e.g., D'Ercole et al. 
2008; Renzini 2008) in which star clusters with masses high 
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enough to retain ejecta in slow winds of stars of the first 
generation gather this material in their central regions where 
second-generation star formation can occur. The observed re- 
lation between the C re i parameter and the initial escape veloc- 
ity of the cluster strongly suggests that this is taking place. 
In the context of this scenario, we recall that the hitherto 
suggested source(s) of the ejecta are FRMS (Decressin et al. 

2007) , massive binary stars (de Mink et al. 2009), and LM- 
AGB stars (e.g., D'Antona & Ventura 2007). Note that the 
ejecta from FRMS and massive binary stars are produced on 
time scales that are significantly shorter than those from IM- 
AGB stars (~ 10-30 Myr versus ~ 50-300 Myr, respec- 
tively; see e.g. Decressin et al. 2007; Ventura & D'Antona 

2008) . I.e., all material created by massive stars that could 
be used to form second generation stars would be available in 
the cluster by the time IM-AGB stars would start their slow 
mass loss. To compare these time scales to the age distribu- 
tions of the eMSTO clusters in our sample, we show the latter 
in Fig. 6, reproduced from Papers I and II. Note that Monte 
Carlo simulations described in Paper II show that two popu- 
lations with ages separated by 100-150 Myr or more would 
result in observable bimodality in our age distributions de- 
rived from the MSTO photometry. Hence, the combination 
of the observed age ranges of 200-500 Myr, the fact that the 
material lost from first-generation stars throughout the cluster 
needs some time to accumulate in the central regions to allow 
second-generation star formation to occur, and the absence 
of clear bimodality in the age distributions seems to suggest 
at face value that FRMS and/or massive binary stars could 
well be significant contributors to the enriched material used 
to produce the second stellar generation. However, IM-AGB 
stars also seem likely significant contributors, since the age 
distributions of the star clusters with the highest escape ve- 
locities in our sample do typically peak near the younger end 
of the observed age range (especially for the more massive 
clusters in our sample). 

Our results also have an impact on an argument made 
against the "in situ" scenario that if indeed the eMSTO phe- 
nomenon reflects an range of ages and the second generation 
of stars is formed from material produced by first-generation 
stars, why are no eMSTOs seen in clusters with ages simi- 
lar to the typical spread in ages within eMSTO clusters (i.e., 
100-500 Myr)? The answer seems to be, at least in part, 
that the escape velocities of clusters in this age range that 
are nearby enough to detect the presence of an eMSTO are 
not high enough to retain the ejecta of the first-generation 
stars. Using the masses and radii derived by the study of 
Mackey & Gilmore (2003) for Magellanic Cloud clusters in 
that age range and applying our methodology described in 
§ 3.3, we find that there is only one such cluster (NGC 1856, 
log Mci/Mq — 4.9, r h ~ 3.2 pc, log age = 8.12) that has a 
v esc ,7 value above 15 kms" 1 , i.e., sufficient to retain gas lost 
by slow stellar winds of the first generation. Unfortunately, 
the currently available photometric data for this cluster are in- 
adequate to properly assess the presence of a range of ages. 
Keller et al. (2011) also recently made an additional, statisti- 
cally based argument suggesting that the absence of eMSTO 
clusters in the age range of 100-300 Myr is not problematic 
to the "in situ" scenario. 

In addition to the "in situ" scenario, part of the gas used 
to form second-generation stars may have been accreted 
from the ambient ISM > 10 8 yr after the first generation 
was formed (see Conroy & Spergel 2011; Pflamm-Altenburg 



& Kroupa 2009). This scenario is attractive in that (;) 
it addresses the concern that the high fraction of second- 
generation stars in GCs seem to require more gas than any 
suggested type of donor stars can provide, and ( ii) gas from 
the ambient ISM readily provides the "dilution" that is likely 
needed to create the smooth Na-O anticorrelations seen in an- 
cient Galactic GCs (e.g., Ventura & D'Antona 2008, 2009). In 
the absence of such dilution, one would expect a more discrete 
distribution of cluster stars in the [Na/Fe] vs. [O/Fe] diagram 
than that observed (see, e.g., Carretta et al. 2009). However, 
the uniformity of massive element abundances (e.g., [Fe/H] 
and [Ca/H]) among stars in Galactic GCs (except for the most 
massive ones like us Cen and M 54) and the razor-sharp RGB 
sequences seen in the CMDs of the more massive LMC star 
clusters in our sample (cf. Paper II) do require that the spread 
in [Fe/H] and [Ca/H] in the ambient ISM was very small when 
the second-generation stars were formed. Future measure- 
ments of [Fe/H] and [Ca/H] variations among RGB stars in 
the LMC field population may be able to address this con- 
cern. A next generation of theoretical simulations similar to 
those used by D'Ercole et al. (2008, 2010) may well be able 
to shed more light on the expected relative contributions of 
the possible donor star types after taking into account their 
respective numbers of stars from the IMF (see, e.g., de Mink 
et al. 2009), their typical evolutionary time scales, the abun- 
dance patterns of their ejecta, and the mass evolution of star 
clusters containing multiple generations of stars. 

5. SUMMARY AND CONCLUSIONS 

We analyzed dynamical properties of 1 1 intermediate-age 
star clusters in the LMC that have high-precision ages and 
metallicities measured from recent HST/ACS photometry. 7 
of these 11 clusters contain MSTO regions that are signif- 
icantly more extended than expected from measurement er- 
rors. 

For the three massive clusters NGC 1751, NGC 1806, and 
NGC 1846, we find that radial distributions of stars in the 
upper (i.e., brighter) half of the eMSTO region are signif- 
icantly more centrally concentrated than stars in the lower 
(i.e., fainter) half of the eMSTO region as well as the (more 
massive) RGB and AGB. Since this cannot be due to dy- 
namical evolution of a SSP, we conclude that the upper and 
lower MSTO regions in those clusters correspond to intrin- 
sically different populations which have undergone different 
amounts of violent relaxation during their collapse. To look 
for an explanation for this finding, we calculate the evolution 
of mass and escape velocity for the clusters in our sample 
from an age of 10 7 yr to the present time, using a combination 
of stellar evolution models and dynamical cluster disruption 
models. The three star clusters mentioned above turn out to 
be the only three clusters in our sample that have the follow- 
ing two properties: (7) escape velocities v esc > 15 kms" 1 at 
an age of 10 7 yr, i.e., large enough to retain material shed by 
slow stellar winds of the first generation, and (ii) half-mass 
relaxation times that are larger than or equal to their age. At 
least for these clusters, we suggest that the main cause of the 
presence of the eMSTO region is a range of stellar ages within 
the cluster. The data seem most consistent with the "in situ" 
scenario in which secondary generations of stars form within 
the cluster itself out of gas ejected by stars of the first stellar 
generation that feature slow stellar winds. Viable sources of 
the enriched material are thought to include fast-rotating mas- 
sive stars, massive binary stars, and intermediate-mass AGB 
stars. Element abundance ratios from high-resolution spec- 
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troscopy of individual cluster stars should be very useful in 
further constraining the nature of the eMSTO regions in mas- 
sive intermediate-age star clusters. In particular, one would 
expect to see correlated variations in light element abundances 
(e.g., N, O, Na) among the stars in the star clusters with rel- 
atively high escape velocities in our sample, likely in a way 
similar to the Na-0 anticorrelation found in ancient GCs. 
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Table 2 

Adopted dynamical parameters of the eMSTO star clusters studied in this paper. 



; logCM c i/M ) I r h ; ( rclax I Vesc ; 

Cluster Current 10 7 yr Current 10 7 yr Current 10 7 yr Current 1() 7 yr C re i 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



Case I: Model Without Initial Mass Segregation. 



4.82 ± 0.09 
5.25 ± 0.09 
5.03 ± 0.09 
5.17 ±0.09 
4.49 ± 0.09 
4.41 ± 0.09 
4.00 ± 0.09 



5.07 ± 0.09 
5.46 ± 0.09 
5.27 ± 0.09 
5.41 ± 0.09 
4.73 ± 0.09 
4.66 ± 0.09 
4.31 ± 0.09 



7.1 ±0.6 
14.7 ±3.2 
8.9 ± 0.6 
8.4 ± 0.4 
8.9 ± 2.2 

7.1 ±0.6 

7.2 ± 1.7 



5.0 ±0.4 
9.9 ± 2.2 
6.3 ± 0.4 
5.9 ±0.3 
6.2 ± 1.5 

5.0 ± 0.4 

5.1 ± 1.2 



1.2 ±0.2 
5.2 ± 1.8 

2.0 ± 0.2 

2.1 ±0.2 

1.2 ±0.4 
0.8 ±0.1 
0.6 ± 0.2 



0.9 ± 0.1 
3.8 ± 1.3 
1.5 ±0.2 
1.5 ±0.1 
0.9 ± 0.3 
0.6 ±0.1 
0.4 ± 0.2 



9.2 ± 0.9 
9.9 ± 1.5 
10.3 ± 1.2 
12.7 ± 1.2 

5.6 ± 0.8 

5.7 ± 0.6 
3.5 ± 0.5 



14.6 ± 1.5 
12.6 ± 2.3 
16.3 ± 1.6 
20.0 ± 1.9 
8.8 ± 1.3 
9.1 ±0.9 
6.0 ± 0.9 



2.52 ± 0.53 
1.33 ±0.07 
2.06 ±0.16 
3.36 ±0.18 
1.55 ± 0.37 
1.45 ±0.62 
1.27 ±0.61 



Case II: Model Involving Initial Mass Segregation. 



4.82 ± 0.09 
5.25 ± 0.09 
5.03 ± 0.09 
5.17 ±0.09 
4.49 ± 0.09 
4.41 ± 0.09 
4.00 ± 0.09 



5.17 ± 0.09 
5.58 ± 0.09 
5.37 ± 0.09 
5.52 ± 0.09 
4.82 ± 0.09 
4.75 ± 0.09 
4.41 ± 0.09 



7.1 ±0.6 
14.7 ±3.2 
8.9 ± 0.6 
8.4 ± 0.4 
8.9 ± 2.2 

7.1 ±0.6 

7.2 ± 1.7 



2.6 ± 0.2 
5.5 ± 1.1 
3.3 ± 0.2 
3.1 ±0.1 
3.3 ±0.8 

2.7 ±0.2 
2.7 ± 0.6 



1.2 ±0.2 
5.2 ± 1.8 

2.0 ± 0.3 

2.1 ±0.3 

1.2 ±0.5 
0.8 ±0.1 
0.6 ± 0.2 



0.4 ± 0.1 
1.6 ±0.6 
0.6 ±0.1 
0.7 ±0.1 
0.4 ± 0.1 
0.3 ±0.1 
0.2 ±0.1 



9.2 ± 0.9 
9.9 ± 1.5 
10.3 ± 1.2 
12.7 ± 1.2 

5.6 ±0.8 

5.7 ± 0.6 
3.5 ± 0.5 



22.6 ± 2.3 

20.2 ± 3.0 

25.3 ±2.5 
31.0 ±3.0 
13.5 ±2.0 
13.5 ±2.0 

9.3 ± 1.3 



2.52 ± 0.53 
1.33 ± 0.07 
2.06 ±0.16 
3.36 ±0.18 
1.55 ± 0.37 
1.45 ±0.62 
1.27 ±0.61 
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Note. — Column (1): Name of star cluster. (2): Logarithm of adopted current cluster mass (in solar masses). (3): Logarithm of 
adopted cluster mass at an age of 10 Myr. (4): Current cluster half-mass radius in pc. (5): Adopted cluster half-mass radius at an age of 
10 Myr. (6): Current cluster half-mass relaxation time in Gyr. (7): Cluster half-mass relaxation time at an age of 10 Myr. (8): Current 
cluster escape velocity at tidal radius in kms~'. (9): Cluster escape velocity at tidal radius at an age of 10 Myr. (10): Value of C re i, 
defined in § 3.3.4. 



